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ABSTRACT 

The IceCube collaboration recently reported the potential detection of two cascade 
neutrino events in the energy range 1-10 PeV. We study the possibility that these PeV 
neutrinos are produced by gamma-ray bursts (GRBs), paying special attention to the 
contribution by untriggered GRBs that elude detection due to their low photon flux. 
Based on the luminosity function, rate distribution with redshift and spectral properties 
of GRBs, we generate, using Monte-Carlo simulation, a GRB sample that reproduce 
the observed fluence distribution of Fermi/GBM GRBs and an accompanying sample 
of untriggered GRBs simultaneously. The neutrino flux of every individual GRBs is 
calculated in the standard internal shock scenario, so that the accumulative flux of 
the whole samples can be obtained. We find that the neutrino flux in PeV energies 
produced by untriggered GRBs is about 2 times higher than that produced by the 
triggered ones. Considering the existing IceCube limit on the neutrino flux of triggered 
GRBs, we find that the total flux of triggered and untriggered GRBs can reach at most 
a level of ~ 1CP 9 GeVcm - 2 s _1 sr _1 , which is insufficient to account for the reported two 
PeV neutrinos. Possible contributions to diffuse neutrinos by low- luminosity GRBs and 
the earliest population of GRBs are also discussed. 

Subject headings: diffuse background : neutrinos — gamma rays: bursts 



1. Introduction 

Detection of high-energy neutrinos would be an important step towards identifying the origin 
of cosmic rays. Interestingly, through analysis designed to search for ultra high-energy (> PeV) and 
high quality shower events over a period of 672.7 live days between 2010 a nd 2012 with IceCube, 



two PeV neutrino-induced cascade events have been found ( Ishihara Il2012l ). The two events were 



observed on August 9, 2011 and January 3, 2012, and each have an energy of approximately 1-10 
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PeV. The observation of two events with an expected background of 0.14 events from the atmo- 
spheric neutrinos corresponds to a significance of 2.36cr. Given this modest statistical significance, 
we cannot at this time be confident that they are of astrophysical origin. 

Gamma-ray bursts (G RBs) have been proposed as a potential s ource of ultra-high energy 
cosmic rays (UHECRs) (e.g. IWaxmanlll995l : IVietrilll995l : iDermerl 12002] ) . and high energy neutrino 
emission has been predicted to be produced in the dissipat i ye fireballs, where cosm i c ray protons 



interact with fireball pho t ons (e.g. |W axman &; Bahcalllll997l ; iDermer fc Atovanl 120061 ; iGuetta et al 



20041 : iMurase et al.ll2006l : IWang fc Daill2009^ . Gamma-ray bursts are an attractive possibility for 



the reported two PeV neutrinos because the predicted neutrino emission in the standard internal 
shock model peaks at PeV energies. Notably, previous searches of GR B neutrinos with IceCube 
via s tacking analysis of many triggered GRBs have yield none detection (|Abbasi et alJl2010l . 



2011 



2012l ). The approach to search for neutrinos from triggered GRBs is through analysis of upgoing 
neutrinos that are correlated in time and direction with known GRBs, which is different from that of 
detecting diffuse neutrinos. The most stringent constraint is given by the combined analysis of the 
IceCube 40- and 59-string (IC40+59) data of upgoing neutrinos fro m 215 triggered GR Bs, which 
leads to an upper limit that is claimed to be below the expected fluxl Abbasi et al. 2012l fl Beside 
the contribution by triggered GRBs for which the stacking analysis has been performed, GRBs that 
do not trigger detectors due to their low flux may also contribute to the diffuse neutrino flux. The 
purpose of this paper is to study whether the sum contribution of triggered GRBs and untriggered 
GRBs can explain the reported two PeV neutrinos under the constraint that the stacked neutrino 
flux of 215 triggered GRBs is below the IC40+59 sensitivity. The untriggered GRBs defined in our 
paper are normal GRBs (i.e. belonging to the normal GRB population) that do not trigger the 
current detectors. They are normal GRBs in the low- luminosity end of the luminosity function (for 
normal GRBs) and/or locating at high redshifts and hence they can elude the detection by current 
GRB detectors (e.g. Fermi/GBM) due to the limited sensitivity. 

Previous estimates of diffuse neutrino flux from GRBs are made through integrating the neu- 
trino flux over the whole luminosity range and the redshift range of GRBs, assuming that the 
their number density follows a given luminosity function an d redshift distribution function (e.g. 
Gupta fc Zhang) 120071 ; Irle et all 12012 ; Icholis & Hooper 112012 ). The estimates show that the total 
flux from all the GRBs (including both the triggered ones and the untriggered ones) reach a leve l 
of ~ 10~ 9 - 10- 8 GeVcm- 2 s- 1 sr- 1 (|Gupta fc Zhangl I2OO7I : \lie et~alll2012l : lcholis fc Hooperlboilj l. 
which depends on the luminosity function and redshift distribution that are used. These estimates 
made simplified assumptions about some properties of GRBs, e.g. assuming the same photon 
spectrum for every GRBs, and do not consider the existing IceCube limit on triggered GRBs. To 
overcome this simplification, we here present a new approach to calculate the diffuse neutrino flux. 



1 This claim is based on the theoretic flux obtained with the formula ir Guetta et al. (2004). Recent calculations 



show , however, that the previous calculation may overestimate the flux of GRB neutrinosfL 
2012l : lHe et ahlboid l. 



2012 



Hummer et al 



- 3 - 



We generate a GRB sample with Monte-Carlo simulation □, requiring that the triggered population 
of GRBs reproduce the observed properties of Fermi/GBM GRBs. A population of dim, untrig- 
gered GRBs will be generated simultaneously in the simulation. Then, the neutrino flux of every 
individual GRBs are calculated, so that the accumulative flux of the whole sample can be obtained. 

The rest part of this paper is arranged as follows. In §2, we describe the properties of GRBs 
and some empirical relations used in our simulation generating the GRB sample. We calculate 
the accumulative neutrino flux from both triggered GRBs and untriggered GRBs in the generated 
sample and compare them with the reported two PeV events in §3. We give our discussion and 
conclusion in §4 and 5. Throughout the paper, we use eV for particle energy, c.g.s units for other 
quantities and denote by Q x the value of the quantity Q in units of 1(F unless specified. 



2. Simulation of a GRB sample 



The distributions of GRBs in luminosity and in redshift are assumed to be independent of 
each other in this work, i.e., n(L,z) = <p(L)p(z), whe re ML) and p(z) are lumino sity function and 
redshift distribution function of GRBs respectively. IWanderman fc Piranl ()2010l . hereafter, W10) 
obtained the luminosity function and cosmic rate by inverting directly the redshift - luminosity 
distribution of observed long Swift GRBs, which are, respectively, given by 



and 
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where a\ = 1.2, a 2 = 2.4, L* = 10 5i - 5 ergs , n\ = 2.1, n 2 = -1.4, = 3.1, and p = l.SGpc^yr -1 
are the best fit parameters. Note that L is the isotropic, peak bolometric luminosity of GRBs in 
energy range 1 keV — 10 MeV. The low- and high- luminosity ends of the luminosity function are 
set as lO^ergs 1 and 10^ 4 ergs 1 respectively, and the maximum redshift is set as z mSLX = 8. We 
generate a sample of GRBs that follow the luminosity distribution and redshift distribution given 
by Eq. [1] and Eq. [2] with Monte-Carlo simulation. The isotropic energy of each GRBs is obtained 
through the relation 

log £ iS o,52 = 1.071og L iso ,52 + (0.66 ± 0.54) (3) 



which is obtained from the data in iGhirlanda et al.l (|2012l ). The photon spectra o f GRBs can be 
described by an empirical function known as the Band function ([Band et al.lll993l ). which can be 



2 A similar app roach of sampling the bursts individually using a Monte-Carlo algorithm was made by 
Baerwald et al. to obtain a more complete sample of GRBs. 
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written as 



N(E) = NO < 
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where Nq is the normalization constant and a, (3 are the photon indices for low- and high-energy 
spectra respectively. Generally, a > — 2 and (3 < —2, and the energy spectrum peaks at -Speak = 
(2 + ol)Eq. The values of a and j3 are generated following the observ ed distribution of BATSE long 
GRBs in the simulation and only bursts with j3 < — 2 are adopted ([Kaneko et ajj 120061 ) . The rest 
frame peak energy of GRBs is found to be correlated with their isotropic ener gy and luminosity , 
which are known as t he Amati relation and the Yonetoku relation, respectively (jAmati et al.l 12002 : 
Yonetoku et al.l 12004). In this work, to generate the peak energies of GR Bs in the simulation, we 



adopt the -E pea k — -E-iso relation recently found by iGhirlanda et al.l (|2012j), i.e. 



log 
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(5) 



26.06 ± 1.14, 



where -E<°g ak is the observer frame peak energy (hereafter the subscript "ob" represents the quantity 
in the observer frame). Since the neutrino emission is produced dominantly by long d uration GRBs, 
we only include bursts with T90 > 2s in our sample. Following iKakuwa et al.1 (120121 ). we determine 
T90 roughly by T90 = (1 + z){E- lso / 'L av e) where L ave = 0.3Li so is the average luminosity of a GRB. 

To compare our simulated sample with observed GRBs by Fermi/GBM, the trigg er threshold in 



our s imulation is set to be the same as that of Fermi/GBM, i.e., 0.74 photons cm 2 s 1 ( Meegan et al 



2009). Taking into account the effect of Earth occultations and the South Atlantic Anomaly (SAA) 
passages, about 65% GRBs above the Fermi/GBM trigger threshold can be observed. To reproduce 
the GBM detection rate of ~ 250 GRBs per year, the value of po is required to be 1.01 Gpc~ 3 yr _1 , 
which is within the la error of the best-fit value obtained by W10. For this value of po, the total 
event rate of GRBs in the luminosity range from 10 50 ergs _1 to 10 54 ergs _1 and redshift range from 
z = 0toz = 8is~ 8000 GRBs per year in the generated sample so the vast majority of them are 
untriggered GRBs. The total number of GRBs occurred from z = to z = 8 per year is obtained 
through N = J 4ttDI(z) ^dD c , where D c (z) = /„* ^ with H(z) = H oy /(l + z)*n M + V A being 
the Hubble constant at redshift z in a flat universe. 

In Fig. 1, we compare the 10-1000keV fluence distributions dP/d(logFy) of GRBs above the 
trigger threshold in our simulated sample with that of 482 Fermi/GBM bursta^l. The K-S test of 
two distributions returns a high accepting probability of 89.5%. The total gamma-ray fluence of 
GRBs in 10-1000keV above the trigger threshold in our simulated sample is 0.0048 erg cm~ 2 when 



3 From GRB 080714A to GRB 100704, data are taken from 



http://heasarc.gsfc.nasa.gov/W3Browse/all/fermigbrst.html 
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the number of bursts is normalized to be 482, which also agree well with the total gamma-ray 
fluence, 0.0044 erg cm -2 , of the observed 482 GBM bursts. 

In Fig. 2 we show the distribution of the cumulative bolometric fluence ^ i i ? 7j j(> jRyih) above 
a certain flux F^^h for the subsample of both triggered and untriggered GRBs. As is shown, the 
total gamma ray fluence of the untriggered GRBs (-Ky jim t = 0.0136 erg cm -2 ) is 2-3 times larger than 
that of the triggered ones (i^ tri = 0.0053 erg cm -2 ). This is mainly due to two aspects. First, the 
number of untriggered GRBs is much larger than that of triggered GRBs, although the luminosity 
of untriggered GRBs is on average much lower than that of triggered GRBs; Second, a fraction 
(about 35%) of bright GRBs occulted by the Earth also contribute to the gamma-ray fluence in 
the untriggered GRBs. This demonstrates that the total contribution of untriggered GRBs can be 
important. 

There is significant var iation in the l umin osity function and the redshift distribut i on of GRBs 
in the literature, see e. g. Liang et al. ( 2007 . hereafter, L07)@ and Guetta &: Piran ( 2007 . here- 



after, G07)H. We find that the GRB sample generated using the above approach with the luminosity 
function and rate redshift distribution given in L07 contain much more triggered GRBs than the 
observed ones, while the G07 luminosity function in combination with the relevant redshift dis- 
tribution gives rise to too few triggered GRBs. In other words, to maintain the rate of triggered 
GRBs of 250/yr by Fermi/GBM, p is required to be ~ 0.2Gpc~ 3 yr _1 and ~ 5Gpc _3 yr _1 for L07 
and G07 respectively, which are beyond the 3a error of the best-fit values (1.12 Gpe yr — 1 and 
0.27 Gpc _3 yr _1 for L07 and G07 correspondingly). So here we adopt the luminosity function and 
redshift distribution of W10 in generating the GRB sample in our simulation. 



4 The luminosity function obtained by L07 is 



dN 
dL-y 



+ 



(6) 



where po = 1.12Gpc yr~ is the local rate of GRBs, and $o is a normalization constant to assure the integral over 
the luminosity function being equal to the local rate po- This luminosity function breaks at L 7 b = 2.25 x 10 ergs -1 , 
with indices ai = 0.65 and Q2 = 2.3 below and above the break. The rate distribution wit h redshift used by L07 in 
obtaining this luminosity function is S(z) = 23e 3 ' 4z /(e 3 ' 4z + 22.0') i|Porciani fc Madaull200ll l. 



5 The luminosity function obtained by G07 has the same form as that of W10, but has different parameter values, 
i.e., po = 0.27Gpc _3 yr _1 , oi = —1.1, a-i = —3.0 and L* = 2.3 x 10 51 ergs _1 . This lumin osity function is obtaine d 
based on the assumption that the rate of GRBs follows the star formation history given by iRowan- Robinson! |l999), 
i.e. 

' 10°' 75z z < 1, 



S(z) 



10" 



z > 1. 



(7) 
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3. Diffuse neutrino emission from our GRB sample 
3.1. Analytical calculation of neutrino flux 



To illustrate the dependence of neutrino flux on the parameters, we first present an analytical 
approximation of calculating the neutrino flux. 



It has been suggest ed that proton s can be accelerated up to e° b ~ 10 20 eV through i nternal 
or external shocks (e.g. IWaxmanl Il995l : IVietril Il995l : IWang et al.1 120081 : Murase et al.l 120081 ) . The 
accelerated protons collide with the fireball photons in the same region and produce charged pions. 
The charged pions decay into four final-state leptons via the processes ir + (ir~) — > v^ l (i? fl )+fi + (fi^) — > 
Vaipn) + e + (e~) + v e ipe) + ^u,{ v u)- In the fireball comoving frame, the fractional energy loss rate 
of a proton with Lorentz factor 7 p via photopion production is 



= 7T2 /. dEa 7T (E)( t (E)E 



dxx 2 n(x), 



(8) 



2 7p JE JE/2j p 

where c is the speed of light, E is the photon energy in the rest frame of the proton and Eq is 
the threshold energy of photonpion production, a w is the cross section of photopion production 
and £ is the inelasticity, x and n(x) are the photon energy and the photon spectrum in the GRB 
comoving frame respectively. When only the A resonance channel for the photopion production is 
considered, the fraction of the proton energy converted into the pion can be approximated by 

-/8-l 
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where T is the bulk Lorentz factor of GRB and R is the internal shock radius which relates with 
the bulk Lorentz factor and variability time 5t through R = 2T 2 c5t. The break energy 6^ in the 
proton spectrum is caused by the break in the photon spectrum. Correspondingly, a break occurs 
in the neutrino spectrum at energy 



ob 
e ub 



7.5 x 10 14 (1 + z 
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(10) 



Additionally, there is another break in the neutrino spectrum caused by the synchrotron cooling of 
charged pions in the magnetic field of the dissipative fireball, which can be written as 



ob 



3.3 x 10 17 (1 + 



(11) 



The final neutrino spectrum (single flavor) of a GRB is approximately a three-section power law 
with two breaks at and respectively, i.e., 

1 f]pE^\ >0 \ 
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where rj p is the ratio between the energy in the accelerated protons and the radiation energy, e Pimax 
and e Pim i n are respectively the maximum and minimum energy of accelerated protons, and -F^boi 
is the bolometric (IkeV-lOMeV) gamma-ray fluence of the GRB. The factor 1/8 comes from the 
fact that half of the produced pions are charged pions and each neutrino takes about 1/4 of the 
pion's energy. 



3.2. Numerical results 

In our numerical calculation, besides the A resonance channel, we also consider the contribution 
from the direc t pion channe l and t he multi pion channel. The cross-section for photopion production 



is taken from iMiicke et al.l (|2000l ) and the branching ratio of each channel is approximated as the 
same in H12. Muon cooling, pion cooling and neutrino oscillation effects are also included in the 
calculation. We calculate the GRB neutrino flux within the standard intern al shock model, i.e. 



we u s e R = 2T c5t, similar to the treatment by the IceCube Collaboration (jAbbasi et al.l 12010. 



2011 



20121 ). We consider two case s of the choice of the b ulk Lorentz factor V, one is the relation 
r = 29.8-EjgQ52 that is obtained by iGhirlanda et al.l ( 2012 ). and another is the benchmark choice, 



i.e. T = 10 2 - 5 for all GRBs. As in H12, we fix 5t oh = 5t(l + z) at 0.01s. We assume the accelerated 
proton spectrum to be a power law with an index of —2. The baryon ratio is taken to be r\ p = 10. 

We calculate the neutrino fluence from each GRB in the generated sample and the aggregated 
neutrino fluence from all the triggered GRBs in the sample is shown in Fig. 3. To compare with the 
IceCube limit from the combined analysis of IC40 and IC59 for 215 GRBs, we have multiplied the 
aggregated fluence by a factor of 215/iVtri, where iV tr i = 250 is the number of triggered GRBs in our 
sample. The aggregated neutrino fluence peaks around IPeV and reach a level of ~ O.lGeVcm -2 
for both choices of T. For comparison, we also show the fluence limit of IC40+59 in Fig. 3 (dashed 
lines). The 90% upper limit on the fluence is obtained by requiring that the upgoing muon neutrino 
event numbej§ observed by IC40+59 should be smaller than Nn m i t = 1.9. The non-detection of 
neutrinos by IC40+59 implies rj p < 20. We note that the aggregated neutrino fluence of 215 GRBs 
in our former work H12 (see Fig. 2 in H12) is a factor of two lower and peaks at higher energy. 
The difference is due to that, in that work the luminosity and bulk Lorentz factor of all the GRBs 
without measured redshift are fixed as 10 52 ergs _1 and 10 2 ' 5 respectively. 

The diffuse neutrino flux can be obtained by multiplying the aggregated neutrino fluence of 
the all-sky GRBs occurred in one year by a factor of l/4-7ryr ~ 2.5 x 10~ 9 s _1 sr _1 . Fig. 4 presents 
the diffuse neutrino flux from all GRBs in our generated sample (solid lines), which consists of 



6 The upgoing muon neutrino number is obtained by summing up the contribution from every individual GRBs in 
the energy range of 0.1-3 PeV. For these upgoing neutrinos, we take 90° — 180° angle-averaged effective area for IC40 
and IC59 in the calculation. The ratio of the number of GRBs calculated with the effective area of IC40 and IC59 is 
set to 117:98, in alignment with the realistic number ratio of GRBs that occurred during IC40 and IC59 operation 
period. 



- 8 - 



the flux contributed by the triggered GRBs (dashed line) and the flux contributed by untriggered 
GRBs (dash-dotted lin e). The top panel corresponds to the case that assumes T = 29.8Ef^ 52 



( Ghirlanda et al.1 12012 ) and the bottom to the case of T = 10 2 ' 5 . In both cases, the total flux 
around PeV energy are ~ 10 _9 GeVcm _2 s _1 sr _1 . The flux produced by the untriggered GRBs is 
about a factor of two higher than that of the triggered ones. 

In Fig. 5, we show the contribution to the diffuse neutrino flux by untriggered GRBs with 
different luminosity separately. As can be seen, GRBs in the intermediate luminosity range (10 51 — 
10 53 ergs _1 ) contribute the largest fraction of the diffuse neutrino flux in both cases of the choice 
of the bulk Lorentz factor. These GRBs do not trigger the detectors either due to occurring 
at relatively high redshifts, and/or due to having high spectral peaks (e.g., E^ ak S> 300 KeV). 
Besides, the occulted bright GRBs also contribute a significant fraction to the total neutrino flux. 



3.3. Comparison with two PeV neutrinos 



The IceCu be Collaboratio n recently reported the detection of two cascade neutrinos in the 

neutrinos. Several 



range 1-10 PeV (|lshihara 1120121 ). which are likely to be electron (or anti-electron 
scenarios have been proposed to explain the origin of these two neutrinos 



Baerwald et al. 



2012; 



Bhattacharva et al. 



2012 



Bar ger et al. 



Roulet et al.ll2012l : ICholis Hooper 1 12012ft . Since 



2012 



the energies of the reported two neutrinos are close to the peak energy of the typical GRB neutrino 
spectrum, we here study the possibility that they originate from diffuse GRB neutrinos. In studying 
the possible GRB origin, one must assure that the stacked neutrino flux of those 215 triggered GRBs 
observed by IceCube do not exceed the IC40+59 limit. The expected number of diffuse neutrinos 
from GRBs in the energy range l-10PeV can be obtained by 



N dm 



10PeV 



IPeV 



dn Ue 



(13) 



diff 



where ' 



diff 



is the number spectrum of diffuse electron neutrinos (in unit of cm 2 s x sr x eV 



and S(e u „) is the Ice Cube exposure of electron neutrinos at ultra-high energies in 2 years (2010- 
2012) (jlshihara 1 120121 ). which includes the contribution by the Glashow resonance. Notably, since 
the analysis approach for diffuse cascade neutrinos is different from that for upgoing neutrinos 
from triggered GRBs, the IceCube effective area is significantly different. Based on the total diffuse 
neutrino flux we obtained above, we get N^s ~ 0.1 for the case of F = 29.8Ef s ^ \ 2 and A^ifF ~ 0.2 
for the case of T = 10 2 ' 5 (rj p = 10 is assumed). Since the non-detection of GRBs by IC40+59 
implies r/ p < 20 (see §3.2), we have A^jiff < 0.2 or 0.4 respectively, which is insufficient to account 
for the reported two neutrinos. There is still a small possibility that these two neutrinos arise from, 
for instance, some very strong GRBs that happen to be occulted by the Earth and hence do not 
trigger the detectors. In addition, strong statistical fluctuation can also possibly lead to detection 
of one or two events. 
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3.4. Effect of variation in luminosity function and redshift distribution 



There is significant variati on in the luminosity f unction <f)(L) and the rate distribution with red- 



shift p(z) in the literature(e.g. iFirmani et al.1 12004 . G07, L07, W10. IZitouni et al. 



2008 



Cao et al 



2011 



The variation may be caused by different approaches in obtaining the luminosity function 
and different spectral properties as well as trigger thresholds being used in the calculation. 

Under the requirement that the total neutrino fluence from the triggered GRBs do not exceed 
the upper limit placed by IceCube, the possibility that GRBs are responsible for the two events 
increases if the fraction of neutrino flux from the untriggered GRBs increases. This requires more 
dim and/or distant GRBs to take place. Therefore, luminosity functions with steep slopes or rate 
distribution with higher rate at high redshift are favorable for producing more dim, untriggered 
GRBs. As the slopes in the luminosity function of L07 below and above the break luminosity 
are shallower than that in W10, adopting this luminosity function will lead to a smaller fraction 
of untriggered GRBs. However, the rate redshift distribution used in L07 suggests higher rate 
of GRBs at high redshift, which, on the contrary, will lead to a higher fraction of untriggered 
GRBs. We estimate that the neutrino flux produced by untriggered GRBs is about 90% of that 
of triggered GRBs for the luminosity function and rate distribution in L07. For the luminosity 
function and rate distribution in G07, this factor increases to 2.5 because of the steep slope in the 
luminosity function of G07 and higher GRB rate at high redshift. Even considering the variation in 
the luminosity function and rate distribution in the literature, the diffuse neutrino flux from GRBs 
may be still insufficient to account for the two neutrinos. 



Discussions 



Other possible sources of dim, untriggered GRBs are low-luminosity GRBs (LLGRBs) and 
Population III (Pop III) GRBs. LLGRBs are GRBs with luminosity < 10 e rgs -1 , which const i tute 



a distinct population fr om the normal GRBs(e.g. ISoderberg et al. 



2006 



2009 



Liang et al. 



2007 



Dai 



Virgili et al.ll2009i ). Although only several cases are observed within ~ 200Mpc so far, which 



is mainly due to their low luminosity, the intrinsic local event rate of LLGR Bs is suggested to be 



about two orders of ma gnitude higher than that of the normal GRBs(e.g. iLiang et al. 



2007 



Dai 



2009 



Virgili et al.ll2009l ). Thus, the large local population coul d compensate for the low energy in 



LLGRBs, and lead to a significant amoun t of n eutrino flux (e.g. iMurase et al.ll2006l ; iGupta Zhang 



2007; Murase et al. 2008 



Liu et al.ll201ll ). In lLiu et al.1 ( 201ll ). we calculated the diffuse neutrino 



flux from LLGRBs for different LLGRB luminosity functions and found that the flux can reach a 
level of ~ 10 _9 GeVcm _2 s _1 sr _1 or higher around PeV energies, depending on the parameters of 
LLGRBs that are used, such as the choice of T. We also found that most LLGRBs are untriggered 
and the energy budget of gamma rays in untriggered LLGRBs can be several times larger than that 
in the triggered ones. Thus more than 0.2-0.4 electron neutrinos in the range 1-10 Pe V could come 



from LLGRBs in principle. However, since the origin of LLGRBs are not clear (e.g. IWang et al 
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20071 ) and the properties of LLGRBs have large uncertainties, we should be cautious about this 
estimate. 



Pop III GRBs are thought to arise from the death of the first stars that locate at high redshift 
(z > 10). It is su ggested that the their l uminosity is typically a few times 10 52 ergs 
duration is > 10 4 s( Meszaros Reesll2010l : iKomissarov &: Barkovll2010l : ISuwa Sz Iokall2011 



1 and their 
so they 



may release energy up to 10 56_57 erg in one burst. At such large distances, their gamma-ray flux 
is below ~ 10~ 8 ergcm _2 s _1 , so they may not be able to trigger current GRB detectors. Thus, 
if Pop III GRBs can produce high energy neut r inos, t hese neutrinos wi ll contribute to the diffuse 
background. It has been shown bv llocco et al.l (|2008l ); iGao et al.l (|201ll ) that the neutrino flux in 



1-10 PeV range can reach ~ 10~ 9 GeVcm _2 s _1 sr _1 , or even ~ 10 -8 GeVcm _2 s _1 sr _1 when some 
preferred parameters are used, e.g., with large progenitor mass and/or high interstellar medium 
density. However, since the flux peaks at ~ 10 17 eV, one would expect even stronger flux at ~ 10 17 
eV, which is inconsistent with the non-detection by IceCube at such en ergies. Moreover, the rate 

is from 
and 



of Pop III GRB s is qu ite uncertain. The Pop III GRB rate adopted in 



Bromm &: Loebl (120061 ) , which obtains a much higher rate than that in Ide Souza et al 



Gao et al 



(2011 



(2011 



Campisi et al.l (|201ll ). So whether Pop III GRBs are responsible for the two neutrinos remains to 
be studied. 

Finally we note that we only consider the neutrino emission under the standard internal shock 
mod el. For other models for the GRB prompt emission, such as the pho tosphere emission m odels 
(e.g. iRees &: Meszaros 1 12005) or magnetically reconne ction models (e.g. Zhang &: Yan 1120111). th e 



efficiency for producing neutrinos / P7 is different (e.g. IWang &: Dai 



2009 



Zhang fc Kumar M2012I ) 



due to different dissipation radii R involved. However, the ratio between the flux from triggered 
GRBs and that from untriggered GRBs is not expected to change. Therefore, the diffuse neutrino 
flux can not be readily increased under the constraint that no neutrinos are detected from from 
triggered GRBs. 



5. Conclusions 

In this paper, we studied the diffuse neutrino emission from GRBs in light of the recent 
report that two PeV cascade events are detected by IceCube. We first generate a GRB sample 
that reproduces the realistic properties of the observed GRBs by Fermi/GBM. An accompanying 
sample of untrigged GRBs is generated simultaneously. Then we calculated the neutrino flux from 
the triggered and the untriggered GRBs in the sample within the standard internal shock scenario. 
The advantage of this approach is that the spectral properties (i.e. the spectral peak and indices) 
of every individual GRBs can be taken into account properly in the calculation. We found that the 
neutrino flux from the untriggered GRBs is about a factor of two larger than that of the triggered 
GRBs. Under the constraint of the non-detection of any upgoing muon neutrinos from 215 triggered 
GRBs in the IC40+59 combined analysis, the total diffuse neutrinos flux (single flavor) can reach 
a level of ~ 10 _9 GeVcm _2 s _1 sr~ 1 at most. Such a flux can only give rise to < 0.2 — 0.4 cascade 
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neutrino events in 1-10 PeV for IceCube 2-years exposure. Low-luminosity GRBs and/or Pop III 
GRBs can also contribute to the diffuse neutrino flux, but large uncertainties in their properties 
prevents us from drawing a firm conclusion. 

We thank Kohta Murase, Jun Kakuwa, Aya Ishihara, Shigeru Yoshida and Ilias Cholis for 
useful discussions. This work is supported by the 973 program under grant 2009CB824800, the 
NSFC under grants 10973008, 11273016 and 11033002, the Excellent Youth Foundation of Jiangsu 
Province and the Fok Ying Tung Education Foundation. 
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Fig. 1. — Comparison of the fluence (10-1000 keV) distribution of GRBs in the generated sample and 
in the observed sample by Fermi/GBM. The solid line and dashed line represent the triggered GRBs 
in our generated sample and the observed sample respectively. The K-S test gives a probability 
of 89.5% for the agreement the two distributions. The inset figure shows the distributions of 
accumulated fluence of these two samples. 
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Fig. 2. — Cumulative bolometric fluence (1-10000 keV) as a function of the flux (50-300KeV) 
of GRBs in the generated sample. The triggered GRBs, untriggered GRBs and the total are 
represented by the dashed line, dot-dashed line and solid line respectively. Since about 35% GRBs 
with flux above the trigger threshold are ascribed to untriggered GRBs due to occultation by the 
Earth and the South Atlantic Anomaly passages, the cumulative fluence of untriggered GRBs is 
non-zero below the Fermi/GBM trigger threshold 0.74phcm _2 s _1 . 
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Fig. 3. — The aggregated neutrino fluence spectrum (single flavor) from triggered GRBs in 
our generated sample. The blue and green solid lines represent the fluence obtained assuming 
r = 29.8E® S Q 1 52 (Ghirlanda et al. 2012) and the benchmark value T = 10 2 ' 5 respectively. The 
corresponding dashed lines represent the 90% upper limits from IC40+59. In the calculation, we 
assumed that the proton spectrum is a power law with index of —2 and ij p = 10. 
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Fig. 4. — The diffuse neutrino flux from GRBs in our generated sample. The dashed line and 
the dash-dotted line represent the contribution by the triggered GRBs and untriggered GRBs 
respectively, and the solid line represnets the total flux. In the top panel T = 29.8Ef s ^ \ 2 is used, 
while in the bottom panel V is fixed at 10 2,5 . We assumed that the proton spectrum is a power law 
with index of —2 and rj p = 10. 
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Fig. 5. — Contribution to the diffuse neutrino flux by untriggered GRBs with different luminosity. 
The dotted lines represent the contribution by the occulted bright GRBs. The dashed lines, dash- 
dotted lines and the long dash-dotted lines represent the contributions from GRBs in different 
luminosity ranges respectively (not including those occulted GRBs). The solid lines represent the 
sum of the four components. ). 



